February 2, 2008 6:26 WSPC/Guidelines stability 



International Journal of Modern Physics D 
© World Scientific Publishing Company 



On stability of simplest nonsingular inflationary cosmological models 
within general relativity and gauge theories of gravity 



G. V. VERESHCHAGIN 

Belorussian State Pedagogical University, Physics Department, 
Sovetskaya str. 18, 220050 Minsk, Belarus 
and 

I.C.R.A. - International Center for Relativistic Astrophysics, 
University of Rome "La Sapienza" , Physics Department, 
P.le A. Moro 5, 00185 Rome, Italy. 
Email: veresh@icra.it 



In this paper we provide approximate analytical analysis of stability of nonsingular infla- 
tionary chaotic- type cosmological models. Initial conditions for nonsingular solutions at 
the bounce correspond to dominance of potential part of the energy density of the scalar 
field over its kinetic part both within general relativity and gauge theories of gravity. 
Moreover, scalar field at the bounce exceeds the planckian value and on expansion stage 
these models correspond to chaotic inflation. Such solutions can be well approximated 
by explicitly solvable model with constant effective potential (cosmological term) and 
massless scalar field during the bounce and on stages of quasi-exponential contraction 
and expansion. Perturbative analysis shows that nonsingular inflationary solutions are 
exponentially unstable during contraction stage. This result is compared with numerical 
calculations. 
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1. Introduction 

In conventional cosmological scenario the Universe has a beginning in time, so called 
Big Bang, associated with initial singularity of Friedmann equations Singular 
state with infinite values of energy density and curvature is unacceptable from 
physical point of view and possibly marks the edge of applicability of General 
Relativity (GR). In spite of some p rogr ess in last decades, attempts to resolve this 
problem within GR had no success 

Correct description of the early Universe when the energy density was very 
high should be given by particle physics 1 At the same time experimental and 
observational data are still insufficient to discriminate between different models of 
elementary particles physics. However in most models scalar fields play an important 
role. One of the main discovery in modern cosmology is a possibility to explain a 
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number of puzzles of Big Bang model by inclusion of short stage of exponential 
expansion of the Universe, so called inflation l*™. During this stage the energy 
density of the Universe is dominated by a scalar field (or fields) and transition to 
usual hot Universe is achieved by creation of relativistic particles at the end of 
inflation El. 

Unfortunately even within inflationary models there is little chance to avoid 
cosmological singularity^. The evolution of the Universe is usually considered from 
the planckian moment when description in terms of classical theory of gravitation 
becomespossible. At contraction stage the singular state is unavoidable in most 
models l^ili 

Nevertheless, some nonsingular inflationary solutions exist Although 
generality of the se solutions i s known to be small they are still of some interest 



in the literature * ' ' ' '. In fact, inflationary initial conditions for the scalar 
field are the same for bouncing inflationary models. It is known that in order to 
have nonsingular cosmological model with the scalar field the potential part of the 
energy density at the bounce must be larger than the kinetic part. These conditions 
lead automatically to subsequent inflation in the case when the scalar field exceeds 
the planckian value at the bounce. Ther efore it was expected that nonsingular 
inflationary models are also quite probable l ^ l -^ l 

Generality of inflation was studied, for instance, inlZISI In Hit was pointed out 
that although infinite number of nonsingular inflationary solutions exist within GR, 
exponential contraction stage of these solutions is unstable. By defining a measure 
on the phase plane where the Hubble parameter equals zero, the probability of 
bouncing inflationary solutions was estimated to be of order unity At 
the same tim e, probability of nonsingular solutions with massive scalar field was 
estimated El to be P ~ 10~ 43 . 

From the other hand, deviations from Riemannian geometry can occur at very 
high energy density. In more general spacetimes GR should be r epla ced by gauge 
theories of gravity (GTG) that are natural generalizations of GRll^l It was shown 
that nonsingular cosmo logi cal solutions within GTG exist in the case of spinning 
matter^! usual matterEj an d matter in the form of scalar field^. Moreover, when 
energy density becomes small enough the deviation from GR becomes negligible 
and at late stages the Universe can be described by Friedmann equations. This 
could allow generalizing standard cosmology and resolving the problem of initial 
singularity. 

In order to analyze the generality of initial conditions for nonsingular cosmo- 
logical models within GR as well as within GTG in this paper perturbations for 
nonsingular solutions are studied. The paper is organized as follows. In the Sec. 2 
we discuss bouncing inflationary models within GR in the case of minimal coupling 
between the scalar field and gravitation. In the Sec. 3 coupling is assumed to be non- 
minimal within the same models. In the Sec. 4 bouncing inflationary models within 
GTG are discussed. In Sec. 5 and 6 perturbations are studied analytically. In Sec. 7 
we present numerical results and comparison with analytical calculations. The last 
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two sections contain discussion and conclusions. 
2. Bounce in inflationary models within GR 

Cosmological Fridmann equations for homogeneous isotropic models in the case of 
single real scalar field ip minimally coupled to gravity read 



k f a\ 2 8n 



:-2 \ ,, I 3M 2 

a Sir 



(1) 



a 3M2 



[<?-V(<p)]. 



where V(ip) is effective potential of the scalar field, a(t) is the scale factor of the 
Universe, M p is planckian mass, k corresponds to the sign of space curvature, and 
the dot denotes time derivative a . 

Conservation law follows from Q as 

» + .(*)# + £-a m 

From 0) it also follows 121 th a t nonsingular solutions exist for the Universe 
with positive curvature [k — +1) for initial conditions at the bouncc b (a(0) = 0, 
a(0) > 0) as 

2 (O) < U(^(0)). (3) 

Solution of (JTJ in WKB approximation ^° and qualitative analysis of dynamical 
system H for the case of massive scalar field (V\ — m 2 (p 2 /2 with m ~ 10 _6 Af p ) 
show, that nonsingular solutions of this kind have practically zero measure. 

The effective equation of state for most solutions approach 'stiff' equation p = p, 
where p is the energy density a nd p is the pressure within GR^ and ultrarelativistic 
equation p — p/3 within GTCl^on contraction stage. 

It is interesting to note that if initial conditions for the scalar field satisfy condi- 
tions ijSl and ip(0) > M p , then nonsingular solutions correspond to chaotic inflation 
l-^on expansion stage. The latter begins right after the transition from contraction 
to expansion. 

Inflationary stage in chaotic scenario lasts until the scalar field reaches planckian 
value ip ~ M p . From that moment scalar field starts to oscillate with damped 
amplitude within GR. In GTG oscillations in scalar field can be accompanied by 
oscillations of the Hubble parameter 1^4 Usually it is assumed that during coherent 
oscillations of the scalar field the energy density of the Universe is transferred into 
the form of ultrarelativistic matter. 

a 7i = c = 1 system is used. 

b It is convenient to measure time from the moment of bounce. 
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3. Nonminimal coupling 

In the previous section the scalar field assumed minimally coup led to gra vitation. 
However, it is necessary for the renormalizability of the theory l ^l | 22 | 23 | that the 
lagrangian contains an explicit coupling term between the scalar field and the scalar 
curvature. The correct form of the lagrangian appears to be 



L 



MpR 



1 



1 



(4) 



16tt 2 ,r 2 
where the semicolon denotes covariant derivative. 

The conformal coupling case corresponds to £ = — 1/6. For arbitrary £ cosmo- 
logical equations arel^" 



k 



a 

- I W, 



(5) 
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The conservation law (0 now generalized as 
+ 3 - £ - 6£ 



<P 



dV 

dip 



= 0. 



(G) 



Nonsingular cosmological model was obtained in based on equations (J5J, 
It was shown, that during expansion this model corresponds to chaotic inflation 
in spite of the fact, that y>(0) < M p . 

Notice, that both on quasi-exponential stages of contraction and expansion as 
well as during transition from contraction to expansion the scalar field changes 
slowly. This corresponds to the following assumption 

ip = ipo = const, (7) 

and as a result 

<p = (p = 0, 



dV 

dtp 



i>\- 
dt 



0. 



(8) 



V = Vq = const. 
In this approximation equations JSJ reduce to 

a 8nV 



3M p 2(l 



(9) 



This is the second Friedmann equation with modified 'newtonian' constant and 
with the same assumptions. Therefore all results obtained for nonsingular models 
described by equations JTJ in approximation 0,© are valid for nonsingular models 
with nonminimal coupling. 
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4. Bounce in inflationary models within GTG 

Cosmological equations of GTG in the case of single real scalar field with the 
effective potential V(ip) in the homogeneous isotropic Universe are ^^^ 



(l - 3f3ip 2 Z- 1 ) - ZpVipZ- 1 
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UP^V'Z- 1 (1 - S^Z- 1 ) ( £ ) - 3/3Z- 1 (V"<p 2 - (V) 2 + ^{V' fZ- 1 ) = (11) 
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where Z(t) = 1 + 0[<?(t) - 4V(<p(t))], V = ™, V" = indefinite 
parameter f3 has dimension of inverse energy density and when p ~ solutions 
of (I10|l . (|lll) deviate from Friedmann solutions. Conservation law follows from 
cosmological equations of GTG. 

Nonsingular inflationary cosmological model was obtained and discu ssed in ^ 
based on these equations. Later, other models were discussed i n EEE31 However, 
they all correspond to chaotic inflation at expansion stage in spite of the presence 
of some features after inflation, such as oscillations of the Hubble parameter. 

Equations l|10|) , Ijllfl can be essentially simplified if we consider exponential con- 
traction or expansion as well as transition from contraction to expansion. During 
these stages assumptions (0),© can be applied so (fTTi f . (fTT|) reduce to similar to 10 
equations 

2 



8ttVo 
2,M 2 ' 

8ttVo 
3Mf 



(12) 



These equations are identical to 0) with ip — and therefore all results obtained 
from the analysis of GR equations in approximation 10,© will be valid in the case 
of GTG in the same approximation. 



5. Zero order equations 

Let us derive equations governing perturbations dynamics during stages of exponen- 
tial expansion and contraction and during transition from contraction to expansion. 
Notice first of all, that our perturbative analysis has nothing to do with analysis of 
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relativistic cosmological perturbations. In effect, we are interested in the following 
problem: what happens with the solution if the initial conditions at the bounce are 
slightly modified. Remind that usual analysis of perturbations deals with spatial 
distribution of perturbations at some moment of time. 

In the following we assume that the scalar field has the effective potential Vy 
with corresponding constraint on the mass of the scalar field, or V2 = A<£>f/4 with 
A = 10~ 12 necessary to obtain small density perturbations after inflation 121 

If we don't take into account the end of inflation and beginning of quasi- 
exponential contraction stage, nonsingular solutions can be studied in approxima- 
tion V ~ Vo = const. However, we prefer to keep the last term in (J2J) as a constant. 
This means we assume that scalar field derivative is nonzero, but the scalar field 
changes very slowly. Reasons for this 'controversial' from the first glance approxi- 
mation will become clear later. 

Further on, we drop term <p 2 in (JJJ, but keep the second term in J5J). With 
parameters values used in inflationary cosmology this causes small error compar- 
ing to numerical solutions. Moreover, the relative difference between exact and 
approximated solutions (the error) reaches the maximum during transition from 
contraction to expansion twice but goes to the asymptotic value during de Sitter 

regimes. Denoting tp = <p and A — i^fr) we arrive to the system of equations 



b=A 2 a, 



4'=-rj - 3-0, 
a 



(13) 



a=b, 



where r\ is a small constant with dimension Mp. It is assumed r\ — m(2Vo) 1 ^ 2 when 
one deals with V\ and 77 = A _1 / 4 (4Vb) 3 / 2 in the case of V2. 
Analytical solution for this system exists and reads 



bo=aoiA smh(At) 

r]sech(At) 3 
^= 12A— 

k a =aoi cosh(At), 



[9 sinh(At) + sinh(3At)] 



(14) 



where aoi is the minimal value of the scale factor at the moment t = and 0(0) = 0. 
The Hubble parameter is H = - and the corresponding solution is 



H = Ataah(At). 



(15) 



Expanding 114f) in series at t — and restricting ourselves to second order we 
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arrive to 



( H Q =A 2 t, 
ifo=-Vot, 



Uo=flo,(l + \A 2 t 2 ). 



(16) 





Fig. 1. Dynamics of cosmological parameters in nonsingular cosmological model with scalar field 
during transition from exponential contraction stage to inflationary stage is shown. The scalar 
field potential is chosen to be Vo = 10 — 3 Mp. Dotted curves correspond to approximate solution 
1161 while firm curves denot exact solution 1141 . 

This solution describes transition from contraction to expansion. Both exact 
and approximate solutions are represented at figd One can see, that both Hubble 
parameter and scalar held derivative ip are approximately linear in time during 
transition. Duration of transition from contraction to expansion t tr can be estimated 
from these expressions since at t — t tr linear regime breaks down. According to l|13fl 
it happens when H = H max = A at 

1/2 



ttr 



2 
A 



3M 2 
2nV 



(17) 



For typical effective potentials used in inflationary cosmology V\ and the 
transition time ranges from 1.2t v to 2 ■ \Q 2 t p for V\ and to 10 5 i p for Va, where t p is 
the planckian time. 



6. Perturbation equations 

Assuming standard perturbative scheme 

(b=b a + b l , 

ip=ip + ipi, 

, a=ao + oi, 



(18) 



C AU quantities here and below are measured in planckian units. 
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where first order corrections b\, %j>\ and a\ are small, and supposing that functions 
bo , t/jQ and ciq are given by (|14|1 the following perturbation equations can be written 



' Vi [9A 2 V + v 2 (^S 6 e - 3S 2 - 1)] ai- 



27A 2 M p 2 



r — [^ Q (2S 2 + l)(6i - AT aai ) - SAfaQn + Aa 0l S l )] , 

, ai=6i, 
where 

Si = sinhL4(> - s)}, C = coshL4(i - s)], 

(20) 

T Q = tanhL4(t - s)], S e = sech[A(t - s)], 

and s denotes the shift in time before the bounce. 

These equations can be integrated only numerically. Even linearized equations 
where terms like tpf or a\ip\ are neglected unfortunately are not solvable analytically 
since their coefficients depend on time d . 

In order to find approximate equations that represent main properties of solu- 
tions of (|19|l we expand time dependent coefficients in l|19|) in series keeping zero 
order quantities only. This allows to reduce system (|19|l to a system with constant 
coefficients. However we keep dependence on s and therefore our approximation is 
valid for every moment of evolution of the model parametrized by s, but during 
short time intervals. 

Notice, that usually quantities like aiV'i and tpf are neglected since they are 
second order with respect to first order quantities denoted with subscript T'. How- 
ever, in the particular problem we deal with, we focus mainly on perturbations in 
the quantity ip and assume that they can be quite large. We will see, that the term 
with V? in © causes asymmetry in perturbations dynamics during contraction. 
The term i/iibi is neglected in the second equation of (|19|l and in denominator of 
this equation a\ is also neglected. 



d Full set of equations 1191 is involved into subsequent numerical analysis since linearized per- 
turbed equations are, in general, equations describing different problem. Numerical calculations 
of differences of two close solutions agree with the 1191 solutions but disagree in general with 
solutions of the linearized system. 



2, 2008 6:26 WSPC/Guidelines stability 



On stability of simplest nonsingular inflationary cosmological models. . . 9 

The new system of equations reads 

' ai + C 2 ipi + C ipf, 

^Cibx + dax + C^u (21) 

di=&i, 



where 



87ra 0l . . . 
C ° = -3M^ C0Sh(As) ' 

Ci = 27 # M3 [9^ 2 Vb + ?7 2 (4sech 6 (A s ) - 3sech 2 (A s ) - 1)] , 

° 2 = ~JlS SecHAs) tanh ( As )( 5 + cosh(2A s )), (22) 

C\ = — -^-sech(As) tanh(As)(l + 2sech 2 (As)), 
Aa i 

C 4 = C 3 Atanh(As), 
C 5 = 3Atanh(As). 

Notice now, that the coefficient C5 is always much greater than C3 and C4. In 
fact, constant A is small since V <C Expanding these coefficients into series 
and keeping second order terms one can check that C5 3> 63,64. This inequality 
is valid for all s except for s = as at that moment C3 = C4 = C5 = 0. Therefore, 
equation for ipi can be decoupled from (|21|l and integration gives immediately 



^1 = 5ij) exp {C 5 t} = Si/j exp \ -7^(671%) 1/2 tanh 



M p 



2s ( 2irV ^ 1/2 ' 



t}, (23) 



where Sij) is initial deviation from solution ipQ. 

This constitutes the main result of the paper and shows that for negative s < 0, 
i.e. for moments of time after the bounce perturbations of scalar held derivative 
converge exponentially with time, so the motion of the scalar held is asymptotically 
stable. On the contrary, on contraction stage with s > perturbations of scalar held 
derivative grow exponentially with time and therefore the motion is unstable. 

Analytical solutions for functions b\ and a\ in the case Co = can also be 
obtained after substitution of ipi from (|23|l but they look quite intricate. The most 
important point concerning these solutions is that both for contraction and expan- 
sion stages they can be represented as 



a\ = hiexp{a\t} + h-iexp^a^t] , b\ = h^expla^t} + h^expla^}, (24) 
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where ft.1,2,3.4 are constants and 0^3 < 0. Constants 012.4 > on contraction stage 
so at late times both 61 and di diverge exponentially, whereas on expansion stage 
«2.4 < 0. One can object to this conclusion as equations H21(l are valid only for short 
time intervals. Comparison to numerical results shows however that perturbations 
diverge indeed and approximate solutions of l(2"T)l correctly illustrate the situation. 

Finally, it is necessary to note that the coefficient Co grows exponentially with 
s and it is positive for s > 0. Therefore, the role of the term ipf becomes crucial far 
from the bounce on contraction stage. 

7. Numerical results 

Our analysis of the dynamics of perturbations results in a strong conclusion, namely 
that perturbations of the Hubble parameter and the scale factor diverge exponenti- 
ally with time on contraction stage. The same happens with scalar field derivative 
perturbations at contraction stage while after Hubble parameter becomes posi- 
tive all functions a\ , 61 and ipi converge exponentially with time and therefore are 
asymptotically stable. 

Analytical expressions for perturbations are, however, solutions of extremely 
simplified perturbative equations, where time dependence in corresponding coeffi- 
cients is neglected. It is interesting to compare our results with numerical solutions 
of full equations We provide this comparison at fig0 Since function 61 is a 

derivative of the scale factor perturbation a\ it gives insight into the dynamics of 
the Hubble parameter perturbations. 

In effect, since 60 an d V'o cross zero at the bounce this will cause problems 
in understanding dynamics of perturbations in 61/60 representation. So functions 
61, cii and i/ji are shown at corresponding figures but not dimensionless quantities 
61/60 etc. 

We have chosen the value of the scalar field potential as Vo = 10 _3 Mp so the 
transition from contraction to expansion takes about t tr ~ 20t p . Perturbations 
are calculated for s — —30 at fig. when transition is already finished so the 
dynamics of perturbations correspond to inflationary stage. It is clear that exact 
solutions mimic the dynamics for the scale factor and its derivative perturbations. 
On the other hand, on inflationary stage perturbations in the scalar field derivative 
converge but more rapidly than numerical calculations give. 

Analytical solution of Ij2f|l fails to reproduce perturbations dynamics during 
transition from contraction to expansion with negative Hubble parameter, as can 
be seen from the fig. 12b, where parameter s = 35 and calculations are made until 
the bounce. Perturbations for the Hubble parameter and the scale factor are in 
good agreement with numerical calculations. On de Sitter contraction stage pertur- 
bations in ipi diverge exponentially as suggested by l|23f) . moreover, divergence rate 
of i|23[) is smaller than one obtained numerically. At the same time, the dynamics of 
scalar field derivative perturbations is incorrect in approximation used, comparing 
to numerical results during transition from contraction to expansion since corres- 
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Fig. 2. Perturbations dynamics at different stages is shown. Dotted curves are exact solutions of 
approximate equations 1211 . whereas firm curves represent numerical solutions of exact perturba- 
tive equations 1191 . From top to bottom figures show perturbations at de Sitter expansion stage 
or inflation (a) de Sitter contraction stage (b) and the same stage but with large if>x (e). 

ponding perturbations diverge while numerical results suggest finite perturbations. 

Interesting phenomenon occurs when perturbations in the scalar field derivative 
are assumed to be large. The term Cotpf in equation l|21Jl plays the crucial role in the 
dynamics of perturbations b\, a\. Equations Ij21|l can be integrated analytically in 
the presence of Co term. This case is shown at figEt- Again analytical and numerical 
results are in good agreement except for ipi function at transition stage. This is in 
fact the most interesting case because perturbations in the Hubble parameter and 
the scale factor diverge exponentially with time but there is an asymmetry so the 
Hubble parameter tends to large negative values. 

In realistic models approximation V ~ const is satisfied during inflation and it 
is also valid during exponential contraction in nonsingular models. While the scalar 
field changes slightly, scalar field derivative grows exponentially and the kinetic part 
of the energy density soon starts to dominate over potential part. This is why the 
effective equation of state on contraction stage tends to p — p in the case of GR 
and p = p/3 within GTG. Therefore, conditions J3J) cannot be satisfied for majority 
of solutions both within GR (with minimal or nonminimal coupling) and GTG at 
the bounce. 

8. Discussion 

Before going to conclusions it is useful to mention briefly the sequence of assump- 
tions used in this paper. 
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• In Friedmann equations we neglect term with scalar field derivative with 
respect to time and assume that the scalar field is approximately constant, 
so the effective potential is flat. 

• In the scalar field equation (conservation law) it is assumed that the deriva- 
tive of the effective potential with respect to the scalar field is nonzero and 
it is a small constant. In zero order this results in good approximation for 
the derivative of the scalar field dynamics comparing to the case of realistic 
effective potentials. 

• Perturbed equations are not linearized and solved numerically. Results 
agree with differences of close solutions with the same initial conditions, as 
for perturbations. 

• Time-dependent coefficients in linearized perturbed equations are expanded 
into series keeping only zero-order terms so coefficients lose dependence on 
time. These approximate equations are integrated analytically. 

• Approximate linearized perturbed equations with quadratic on scalar field 
derivative term are integrated analytically in the case of large perturbations 
of the scalar field derivative. The result is in good agreement with numerical 
calculations. 

The results suggest that nonsingular solutions with chaotic inflation are expo- 
nentially unstable on contraction stage and therefore are of little interest. In effect, 
our conclusion supports one obtained long time ago in^Elon the base of different 
approaches. 

Behavior of flat inflationary models with massive scalar field within GTG is 
analysed in where it is pointed out that bouncing models have zero measure like 
in the case of GR. Results of this paper support conclusion about instability of such 
models. 

9. Conclusions 

Perturbative analysis of nonsingular inflationary cosmological equations of GR (in 
the case of minimal as well as nonminimal coupling) and GTG is given. In ap- 
proximation of slowly varying scalar field cosmological equations can be reduced 
to the case of massless scalar field and cosmological constant. Approximate equati- 
ons for perturbations with time-independent coefficients are solved analytically and 
compared to numerical solutions of exact perturbative equations. 

It is shown, that perturbations of the Hubble parameter, scale factor and the 
scalar field derivative are asymptotically (exponentially) stable during inflation. 
On the contrary, on exponential contraction stage all perturbations are unstable 
and diverge exponentially. For large perturbations of the scalar field derivative 
and far from the bounce perturbations in the Hubble parameter tend to large 
negative values while for small perturbations of the scalar field derivative the Hubble 
parameter tends equally to large positive and negative values. 
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